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ABSTRACT 


The  Development  of  an  Integrating  Cavity  Absorption  Meter  to 
Measure  Optical  Absorption  of  Pure  Waters  and  Suspended 
Particulates.  (May  1990) 

Robin  Merl  Pope.  B.S..  United  States  Military  Academy 
Chair  of  Advisory  Committee:  Dr.  Edward  S.  Fry 


The  Integrating  Cavity  Absorption  Meter  (ICAM)  provides  an 
instrument  capable  of  measuring  optical  absorption  independent  of 
scattering  effects.  The  measurement  of  optical  absorption  has 
always  been  complicated  by  scattering  effects.  The  most  common 
and  perhaps  simplest  method  of  measuring  absorption  is  based  on 
transmission  of  light  through  the  sample.  Scatterers  produce 
systematic  errors  which  prevent  simple  transmission  type 
absorption  measurements.  The  ICAM  is,  in  principal,  rigorously 
independent  of  scattering  effects.  In  this^pc^tf  the  theory  of  the 
device  is  reviewed,  and  five  developmental  models  of  the  device  are 
described.  Measurement  results  are  presented  fon  absorption 
results  in  the  presence  of  varying  concentrations  of  scatterers; 
absorption  spectra  of  various  pure  water  samples  compared  to 


accepted  values;  and  absorption  spectra  for  various  phytoplankton 

cultures.  Absorption  values  as  low  as  .(X)l/meteni^^^ve  been 
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CHAPTERI 

INTRODUCTION 


Accurate  measurement  of  the  optical  absorption  coefficient  for 
both  pure  and  natural  waters  has  been  complicated  by  the  effects  of 
scattering  particulates.  E.  S.  Fry  has  proposed  a  method  of 
measuring  absorption  Axdiich  eliminates  the  systematic  effects  of 
scattering.  This  device,  hereafter  referred  to  as  the  Integrating  Cavity 
Absorption  Meter  or  ICAM,  offers  the  advantages  of  being  simple, 
rugged,  and  capable  of  in  situ  measurements.  The  theoretical  basis  is 
soimd  and  well  documented.  ^  Five  versions  of  the  ICAM  have  been 
constructed,  calibrated,  and  tested.  Optical  absorption  spectra  of 
pure  water  throughout  the  visible  compare  favorably  with  accepted 
values^. 


Journal  model  is  Applied  Optics 
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CHAPTER  n 

REVIEW  OP  UTERATURB 


Methods  of  measuring  optical  absorption  include:  (1)  carefully 
preparing  a  pure  sample,  usualfy  repeated  distillation  in  quartz, 
and  then  passing  a  beam  of  light  throu^  a  long  tube  filled  with  the 
sample^  (2)  measuring  the  attenuation  of  a  beam  of  light  passed 
through  a  long  tube  of  water,  then  meastirlng  the  scattering 
coefficient  of  the  water  and  mathematical^  determining  the 
absorption  coefficient^  (3)  Using  a  long  tube  with  a  highly  reflecting 
lining  to  minimize  scattering  losses^'  ®  (4)  Using  a  point  source  of 
radiation  and  measuring  the  attenuation  at  a  given  radial  distance  and 
over  a  known  detector  area,  assuming  scattering  in  the  sample  is 
isotropic^  (5)  placing  the  sample  at  the  entrance  window  of  an 
integrating  cavity®  (6)thermal  lei^e  or  photothermal  techniques 
based  on  the  temperature  dependance  of  the  refractive  index® 
(7)adiabatic  laser  calorimetry (8)  split-pulsed  laser  methods  to 
measure  extinction  and  scattering,  an  advanced  version  of  (2)  above. 
(9)  Optoacoustic  techniques  to  measure  sound  waves  produced  by 
thermal  expansion  in  a  water  column  subjected  to  an  incident  light 
beam:  the  expansion  depends  on  the  radiation  absorbed  by  the 
water.** 

The  concept  of  using  an  integrating  caviiy  to  aid  in  measuring 
absorption  has  been  used  by  various  researchers  within  the  last  two 
decades.  Elterman  described  the  use  of  an  integrating  sphere  to 
measure  optical  absorption.*®  Integrating  cavities  have  been  used  to 
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measure  absorption  of  particles  placed  on  the  cavity  floor,  in 
suspension,  or  concentrated  on  fllters.^^.i®,^®  Using  a  similar 
principal  (isotropic  radiation)  A.D.  Clark  has  developed  a  variation  of 
the  integrating  plate  method  (also  known  as  opal  glass  plate  method) 
which  he  refers  to  as  the  integrating  sandwich  to  enhance  absorption 
effects  of  particulate  matter  deposited  on  glass  filters. 
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CHAPTERm 

THEORY 


The  ICAM  uses  a  ^tem  of  integrating  cavities  to  produce 
isotropic  illumination  of  the  sample,  i  To  understand  why  the  ICAM  is 
independent  of  scattering,  let  us  first  consider  the  basic  principle  of  a 
sin^e  integrating  cavity. 

The  inside  stuface  of  a  typical  integrating  cavity  has  a  high  diffuse 
reflectance  (Lambertian).  A  beam  of  radiation  introduced  into  the 
cavity  and  incident  upon  this  Lambertian  surface  is  diffusely  reflected 
into  all  directions.  After  numerous  reflections,  the  interior  of  the 
cavity  is  filled  with  a  very  uniform,  isotropic  radiance  distribution. 

The  radiance  level  is  directly  proportional  to  the  total  optical  power 
introduced  into  the  cavity.  Any  object  placed  inside  the  cavity  will  be 
uniformly  illuminated  from  all  directions.  If  nonabsorbing  scatterers 
are  placed  in  the  cavity,  the  net  Irradiance  incident  on  the  wall  of  the 
cavity  wUl  remain  unchanged.  The  cavity  will  contain  the  same 
amount  of  energy  both  before  and  after  the  addition  of  the  scatterers; 
there  is  no  loss  of  radiant  energy  due  to  scattering.  Even  if  particles 
in  the  sample  scatter  radiation  through  some  preferred  angle,  the 
Lambertian  surface  of  the  cavity  will  uniformly  distribute  the  radiation 
so  that  the  irradiance  striking  the  wall  of  the  cavity  will  be  uniform.  If 
the  scatterers  in  the  cavity  have  some  characteristic  absorption,  the 
irradiance  incident  on  the  wall  of  the  cavity  will  vary  inversely  with 
this  absorption.  By  measuring  this  irradiance  we  can  determine  the 
absorption  Independent  of  scattering  within  the  cavity. 


In  designing  a  ^stem  to  produce  isotropic  illumination  in  an 
integrating  cavity  we  want  to  avoid  any  bright  spot  at  the  point  where 
li^t  is  introduced  into  the  cavity  and  we  want  to  avoid  placing  any 
object  in  the  direct  path  of  this  input  beam.  We  achieve  this  by  using 
a  ^tem  of  integrating  cavities  and  by  careful  selection  of  materials 
for  construction  of  these  cavities.  Our  system  of  integrating  cavities 
consists  of  one  cavity  A  completely  contained  within  another  B  as 
shown  in  Figure  1. 


Flgnre  1.  Cross  section  of  the  integrating  cavity  ^tem.  Our 
device  has  cylindrical  symmetry  about  an  axis  perpendicular  to  the 
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The  cells  A*  and  B*  which  constitute  the  walls  of  these  cavities  are 
made  of  Spectralon^**,  a  material  with  very  high  diffuse  reflectivity 
throu^out  the  visible  and  which  is  easily  machineable.  Cavity’  B  is  an 
air  space.  We  place  our  sample  in  cavity  A.  When  light  enters  cavity 
B  it  will  bounce  between  the  two  cell  walls  filling  cavity  B  with 
isotropic  radiation.  Irradiance  incident  on  the  wall  of  cell  B*  from 
within  B  is  designated  F2. 

The  Spectralon  wall.  A',  allows  some  portion  of  the  incident 
radiation  to  pass  through  it  and  into  cavity  A.  The  portion  of  incident 
radiation  that  is  transmitted  is  proportional  to  (1-p),  where  p  is  the 
reflectivity  of  the  wall.  At  each  reflection,  only  a  very  small 
percentage  of  the  incident  light  enteia  cavity  A.  but  the  radiation 
bounces  around  many  times  in  cavity  B  and  effectively  "tries'*  many 
times  so  that  the  radiance  in  A  builds  up  to  an  appreciable  value. 

Since  the  wall  A*  is  sunoimded  by  a  uniform  radiance,  light  enters 
cell  A  and  the  sample  from  every  direction  thus  avoiding  any  bright 
spot  from  a  beam  entering  the  cavity.  This  unique  dual  cavity  design 
is  one  of  the  distinguishing  features  of  this  instrument. 

Another  Important  feature  is  the  extremety  long  path  length 
achieved  inside  the  integrating  cavity.  In  order  to  measure  very  small 
absorptions  one  must  use  a  long  path  length  through  the  sample.  The 
high  reflectivity  of  the  integrating  cavity  wall  produces  a  path  length 
that  is  many  times  the  cavity  diameter.  For  an  empty  cavity  with  a 
wall  reflectivity  of  .992  (the  reflectivity  of  spectralon  at  550  nm)  light 
will  reflect  125  times  before  reflection  losses  decrease  the  intensity 
by  a  factor  of  e*^.  For  a  cavity  diameter  of  3  cm  this  would  produce  a 
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path  length  of  3.75  meters  and  a  diameter  of  10  cm  would  produce  a 
path  length  of  12.5  meters.  It  is  easy  to  see  that  the  high  reflectivity 
of  the  cavity  walls  gives  us  a  path  length  that  will  be  sensitive  to  very 
small  absorptions. 

Assuming  the  radiance  distribution  within  each  integrating  cavity 
to  be  Isotropic  and  homogeneous,  we  denote  the  Irradlance  on  the 
inner  wall  of  cell  A*  from  within  A  as  Fq.  and  firom  within  A*  as  F^.  If 
the  absorption  of  the  Spectralon  and  the  sample  was  precisely  zero, 
the  equilibrium  condition  would  be  F2-F|sFq.  In  practice  this  is  not 
the  case,  but  we  can  use  conservation  of  energy  to  determine  their 
relationships.  At  the  boundary  between  cell  A*  and  the  sample  in 
cavity  A  .  it  is  clear  that  the  energy  entering  the  sample  equals  the 
energy  leaving  the  sample  plus  energy  absorbed  by  the  sample^. 


Energy 

Energy 

Lost 

Lost 

Energy 

Throu^ 

Through 

Energy 

Supplied 

Wall 

Detector 

Absorbed 

FiAo(l-p) 

-  FoAoCl'p) 

♦  FoAd  ♦ 

4aFoVo, 

where  p  denotes  the  reflectivity  at  the  boundary  (wavelength 
dependent).  Aq  is  the  area  of  the  boundary  surface  of  cavity  A.  Aq  is 
the  area  of  the  fiber  optic  detector  for  measiarlng  Fq.  Vq  denotes  the 
volume  of  the  sample,  and  a  denotes  the  absorption  coefficient  of  the 
sample.  Solving  this  equation  for  the  absorption  coefficient  gives, 

a-^[^Ao(l-p)-{A.(l-p)*A„}].  {2) 
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Now  the  measured  signals  vdiich  are  proportional  to  and  Fq  are 
designated  and  Sq.  respectively.  The  ratio  F2/F0  can  then  be 
replaced  S^/Sq.  If  we  now  also  note  that  for  a  fixed  wavelength  and 
volume,  the  parameters  p.  Aq,  Aq.  and  Vq  are  all  constants,  then  we 
can  write. 

o-K,|l-K,.  (3) 

or.  e3q>licitly  showing  the  dependence  on  Vq. 

“V.  =  K,|l-K.,  (4) 

where,  K^.  K2.  K3,  and  K4  are  wavelength  dependent  calibration 
constants  vdiich  must  be  experimentally  determined. 

Methods  of  determining  the  constants  will  be  discussed  in 
chapter  V. 


CHAPTER  IV 


EXPBROIENTAL  APPARATUS 


Figure  2  shows  a  schematic  representation  of  the  ICAM  system. 
We  use  an  Oriel  model  66005  arc  lamp  system  with  a  150  watt  Xenon 
lamp  as  our  light  source.  A  set  of  thlrty^one  lOnm  bandwidth 
Interference  filters  (1  filter  every  10  nm  from  400  to  700nm} 
determine  the  wavelength.  Our  light  source  is  chopped  by  an  EG&G 
Princeton  Applied  Research  Model  197  li^t  chopper.  Six  plastic 
optical  fibers.  1.0  mm  diameter,  carry  the  U^t  to  the  integrating 
cavities;  three  more  fibers  sample  the  irradiances  Fq.  F2.  and  F2. 
These  latter  fibers  are  coupled  to  Burle  4840  photomultiplier  tubes 
which  send  current  signals  to  a  lock  in  amplifier  referenced  to  the 
light  chopper  frequency.  A  Keithley  199  digital  volt  meter  converts 
the  analog  voltagesCSQ.Si.  and  S2)  from  the  amplifier  into  digital 
signals  and  sends  them  via  a  IEEE-488  bus  to  a  Macintosh  11 
minicomputer.  Using  LabVlew  software  a  virtual  instrument  then 
displays  the  value  of  the  absorption  coefficient 

Five  variations  of  a  sfystem  of  integrating  cavities  have  been 
produced  and  evaluated.  These  will  be  referenced  as  models  I-V. 

This  series  of  models  provides  a  thorough  understanding  of  the 
capabilities  and  limitations  of  the  ICAM  ^tem.  Therefore  it  is 
important  to  describe  each  model  as  well  as  its  associated 
experimental  results. 
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Interference  Light 

Filters  Chopper 


Fl^iire2. 


Schematic  of  the  ICAM  system. 
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A.  Model  I 

Figures  3  and  4  show  cross  sections  of  the  ICAM  integrating 
cavities  for  model  I.  This  ^tem  is  rotationally  symmetric  about  the 
Z-axis.  The  center  cavity  is  made  of  Spectralon  with  a  wall  thickness 
of  6  mm,  outer  diameter  of  5  cm,  and  length  of  10  cm.  The  inner 
cavity  holds  an  83  ml.  sample.  The  outer  cavity  is  made  of  stainless 
steel  with  the  inner  surface  coated  with  Duraflect™.  The  inside 
dimensions  of  the  outer  cavity  provide  a  1  cm  space  between  the 
outer  and  inner  cavities.  Eight  spectralon  dowels  support  the  inner 
cavity  and  position  it  in  the  center  of  the  outer  cavity.  Six  of  these  are 
positioned  in  the  X-Y  plane,  and  2  are  positioned  along  the  Z-axis  at 
each  end  of  the  inner  cavity.  Samples  can  be  exchanged  through  a  3 
mm  diameter  hole  in  the  top  of  the  cavity  which  runs  through  the 
upper  supporting  dowel.  This  exchange  is  accomplished  using  a  long 
nosed  pipet  capable  of  reaching  the  bottom  of  the  inner  cavity.  Six 
optical  fibers  suppty  light  to  the  outer  cavity.  Three  of  these  are 
positioned  at  each  end  of  the  system  and  th^  enter  the  cavity  at  an 
angle  of  60  degrees  firom  normal  to  the  cavity  wall  in  the  Z  direction. 
The  optical  fibers  used  to  measure  Fq,  F^,  and  F2  are  positioned  at  the 
midpoint  of  the  cylinder  axis  and  are  symmetricalty  spaced  every  120 
degrees  in  the  x-y  plane  (figure  4).  The  F2  detector  fiber  terminates 
at  the  ixmer  surface  of  B*.  The  Fj  detector  fiber  terminates  at  the 
midpoint  of  the  inner  wall.  A'.  The  Fq  detector  fiber  terminates  at 
the  inner  surface  of  A'.  The  Fq  and  Fj  fibers  are  shielded  firom  the  F2 
Irradlance  by  vin^  tubing  which  is  then  coated  with  a  vdiite  silicone 


^ue.  All  9  fibers  are  shielded  black  vin^  tubing  fit>m  the  exterior 
of  B*  to  their  termination/ongln. 


2 


Input  Fibers 


Figaro  3. 


Model  I  cross  section  parallel  to  the  Z-axis. 


13 


y 

- ►  X 

Flciire4. 

B.  Model  n 

Figures  5  and  6  show  cross  sections  of  the  ICAM  integrating 
cavities  for  model  n  parallel  and  transverse  to  the  Z-axis  respectively. 
Both  cavities  are  15  cm  long  Spectralon  cylinders.  The  outer  cylinder 
is  10  cm  in  outer  diameter  with  a  1.27  cm  wall  thickness.  The  inner 
cylinder  is  5  cm  in  outer  diameter  and  6  mm  wall  thickness. 
Spectralon  end  caps  seal  the  2  cavities  and  hold  the  inner  cell 
coaxially  with  respect  to  the  outer  cell  (figure  5).  Six  input  fibers 
firom  the  arc  lamp  pass  through  the  end  caps  and  terminate  in  the 
outer  cavity.  The  three  fibers  entering  each  end  cap  are  120  degrees 
apart  (m  a  6.4  cm  diameter  circle  and  those  at  one  end  are  offset  by 
60  degrees  with  respect  to  those  at  the  other.  The  inner  cavity  is 
lined  with  a  quartz  bulb  vdiich  tapers  into  quartz  capillary  tubing  at 
each  end.  This  permits  the  exchange  of  samples  without  opening  the 
cavities  and  avoids  direct  contact  between  the  Spectralon  and  the 
sample.  This  quartz  bulb  holds  1 18  ml. 


Model  I  cross  section  transverse  to  the  Z-axis. 
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The  optical  fibers  used  to  measure  Fq,  F^,  and  F2  are  positioned  at 
the  midpoint  of  the  (ylinder  axis  and  are  symmetrically  spaced  every 
120  degrees  in  the  x-y  plane  (figure  6).  The  F2  detector  fiber 
terminates  at  the  inner  surface  of  B*.  The  Fj  detector  fiber 
terminates  at  the  midpoint  of  the  inner  wall.  A'.  This  fiber  is  isolated 
from  the  radiances  in  cavity  B  and  in  the  wall  B*.  by  wrapping  the  fiber 
with  aluminum  foil  and  then  inserting  it  into  a  3/8”  OD  Spectralon 
tube.  To  measure  the  irradiance  from  within  the  sample  cavity.  A.  a 
tube  of  Spectralon  with  1/4”  ID  and  3/8”  OD  is  inserted  into  the  side 
of  the  cavity.  A  detector  fiber  at  the  end  of  this  tube  gives  a  signal 
proportional  to  the  required  irradiance.  To  prevent  the  radiances  in 
A*.  B.  and  B*  from  penetrating  this  tube  we  shield  it  with  aluminum 
foil  and  then  surround  it  with  another  Spectralon  tube  3/8”ID  by  3/4” 
OD.  The  outer  Spectralon  tube  prevents  the  foil  shielding  from 
appredabty  perturbing  these  radiances.  This  ^tem  is  placed  in  a 
light-tight  enclosure  and  the  photomultiplier  tubes  are  mounted 
directfy  on  the  sides  of  the  enclosure.  All  9  fibers  are  shielded 
between  the  outside  wall  of  B'  and  their  termination/origin  points  as 
in  model  I. 

C.  Model  m 

Model  in  is  the  same  design  as  model  n  with  the  exception  that 
the  diameter  of  the  quartz  bulb  is  2.54  cm  leaving  a  air  gap  of 
approximatety  .6  cm  between  it  and  the  inner  wall  of  the  inner 
cylinder.  This  quartz  bulb  holds  75  ml.  Figure  7. 
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D.  Model  IV 

Model  IV  is  identical  to  model  III  with  the  exception  of  a  1.27  cm 
long  band  (or  ring)  made  of  spectralon  positioned  around  the 
midpoint  of  the  quartz  bulb.  This  band  is  approximately  3  mm  thick 
at  its  midpoint  and  tapers  toward  each  end.  Figure  8. 


Figure?. 


Model  ni  cross  section  parallel  to  the  Z-axis. 
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Figure  8.  Model  IV  cross  section  parallel  to  the  Z-axis. 

E.  Model  V 

Model  V  is  shown  in  figures  9  and  10.  Ihis  model  was  designed 
to  hold  a  sample  volume  of  615  ml.  in  a  quartz  glass  bulb  similar  to 
the  one  described  in  ^tem  n.  The  inner  and  outer  cavities  are 
made  of  spectralon.  The  outer  cavity  is  1.27  cm  thick,  25.4  cm  long, 
and  has  a  hexagonal  cross  section  vhen  viewed  along  the  z>axis.  The 
inner  cavity  is  a  cylinder  6  mm  thick,  16  cm  long  and  10  cm  in 
diameter.  The  inner  cavity  is  held  at  the  midpoint  of  the  outer  cavity 
by  spectralon  tubes  of  2.54  cm  OD,  8  mm  ID  and  5  cm  long 
positioned  at  each  end  of  the  ^tem.  Three  spectralon  dowels  at 
each  end  of  the  cavity  8  mm  in  diameter  by  5  cm  long  help  stabilize 
the  inner  cavity.  Input  and  detector  fibers  are  positioned  and 
shielded  as  described  in  model  II. 
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CHAPTER  V 

CALIBRATION 

Two  methods  of  determining  the  calibration  constants  have  been 
used.  We  can  prepare  several  samples  containing  varying 
concentrations  of  an  absorbing  pigment  or  dye.  A  spectrophotometer 
Is  then  used  to  determine  the  absorption  of  the  dye.  The  dye  selected 
for  this  procedure  was  Irgalan  Black^^  (Clba-Gelgy  Corp.)  due  to  its 
falrty  constant  absorption  throughout  the  spectrum  of  Interest  Figure 
1 1  shows  the  absorption  resiolts  for  various  samples  measured  In  a 
Cary  model  219  spectrophotometer.  The  spectrophotometer 
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Flgnze  11.  Absorption  of  Irgalan  Black  samples  measured  in  a 
Caiy  219  spectrophotometer. 
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measures  the  absorption  of  the  dye  relative  to  a  sample  containing 
pure  water  with  no  dye.  therefore  we  do  not  know  the  total 
absorption  of  these  samples  fTotal  absorption  »  absorption  of  dye  + 
absorption  of  water)  but  we  do  know  the  change  in  absorption 
between  each  of  the  samples.  When  placed  in  the  ICAM  each  of  these 
samples  will  produce  a  unique  value  for  the  ratio  (Sj/Sq)  By  plotting 
the  change  in  absorption  vs  (Si/Sq)  and  using  the  method  of  least 
squares  fit,  we  can  determine  the  value  for  Kj  of  equation  (3)  as  the 
slope  of  the  line  fitted  to  these  points.  Figure  12  shows  the  results  of 
this  method  at  a  wavelength  of  630  nm.  Note  that  as  the  dye 


Sl/SO 

Flgnre  12.  Absorption  vs  the  ratio  (Si/Sq)  at  630  nm. 
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concentration  increases  the  Sq  signal  ^will  decrease  linearty.  At 
sufflcientty  high  dye  concentrations  the  absorption  in  the  sample  is  so 
large  that  the  radiation  distribution  becomes  anisotropic  and  the 
decrease  becomes  non-linear.  Specifically,  the  energy  density 
decreases  towards  the  center  of  the  sample  (photons  have  a  high 
probability  of  absorption  before  reaching  the  center).  This  results  in 
absorption  becoming  a  nonlinear  function  of  (Sx/Sq)  as  noted  at  the 
extreme  right  data  point  in  figure  12.  This  Information  also 
establishes  the  dynamic  range  for  the  ICAM  as  the  absorption  range 
that  is  linear  with  respect  to  (Si/Sq).  After  values  for  have  been 
determined,  we  must  determine  the  intercept  value  (K2)  for  equation 
(3).  To  fix  the  intercept  we  must  have  one  sample  where  the  total 
absorption  is  known.  We  measure  the  value  of  (Si/Sq)  with  the  cavity 
empty  and  set  the  absorption  equal  to  zero.  We  can  then  solve  equation 
(3)  for  K2.  The  procedure  Just  described  assumes  the  optical 
properties  of  the  ICAM  do  not  change  appreciably  when  the  inde;<c  of 
refimrtlon  within  the  sample  volume  changes. 

By  determining  the  intercept  with  the  cavity  empty,  the  meter  will 
measure  the  total  absorption  of  axty  solution/suspension  placed  in  cell 
A.  The  absorption  of  suspended  particulates  can  be  determined  by 
measuring  the  sample  with  the  particles  in  suspension,  filtering  the 
solution  to  remove  the  particles,  and  then  measuring  the  filtrate. 
Subtracting  the  filtrate  absorption  from  the  suspension  absorption  will 
give  the  absorption  of  the  particulate  matter  only. 

It  should  be  noted  that  for  measuring  the  absorption  of  suspended 
particulate  matter  or  dissolved  substances,  where  the  absorption  of  the 
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filtrate  must  be  subtracted  firom  the  total  absorption  within  the  cavity, 
the  calibration  procedure  can  be  simplified.  We  can  set  the  absorption 
of  the  filtrate  to  zero  (provided  the  calibration  samples  use  the  same 
medium).  This  allows  both  K^.  and  to  be  determined  from  the  best 
fit  of  absorption  vs  S^/Sq  for  the  calibration  samples  where  one  sample 
contains  pure  medium  (no  added  pigment).  The  absorption  measured 
in  the  ICAM  under  these  constraints  will  be  onfy  the  absorption  of  the 
suspended  particulate  matter.  This  is  equivalent  to  determining  a 
baseline  vdien  using  a  spectrophotometer.  We  then  measure 
absorptions  compared  to  the  baseline.  The  ICAM  offers  the  ability  to 
measure  absorption  of  suspended  particulates  as  well  as  the  absolute 
absorption  of  the  medium,  both  of  vdiich  are  not  possible  in  a 
spectrophotometer. 

In  this  work,  calibration  was  also  accomplished  using  only  one 
solution  whose  absorption  was  accurate^  determined  with  a 
spectrophotometer.  Figure  13  shows  the  measured  ratio.  S^/Sq,  as  a 
function  of  the  volume  of  this  solution  in  the  sample  cavity.  A  least 
squares  procedure  was  used  to  fit  this  data  to  a  straight  line  and  the 
calibration  constants  were  determined  from  the  values  of  the  slope  and 
intercept  using  equations  (3)  and  (4). 

This  latter  calibration  procedure  is  not  generally  used 
because  it  introduces  an  air-liquid  interface  within  the  sample 
volume  and  the  position  of  this  interface  changes  every  time  the 
volume  is  changed.  The  preferred  method  (setting  the 
absorption  equal  to  zero  with  the  cavity  empty  and  using  several 


samples  of  varying  absorption)  always  fills  the  sample  volume 
with  one  homogeneous  sample. 


Flgnre  13.  (S^/Sq)  vs  volume  for  fixed  sample  absorption. 
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CHAPTER  VI 

BIEASUREBSENT  RESULTS 

A.  Independence  of  Scattering 

Of  greatest  Importance  is  the  effect  on  the  observed  absorption 
produced  by  the  addition  of  nonabsorbing  scattering  particulates  to  the 
sample.  Magnesium  Oxide  was  used  for  some  of  the  eaiiy  studies 
conducted  with  model  I.  Ihe  MgO  was  abandoned  due  to  its  solubility 
in  water  and  lack  of  accurate  absorption  data  for  MgO  in  the  visible. 
Quartz  powder  was  selected  as  a  better  alternative  for  the  scattering 
particles.  The  particles  are  irregular  and  of  the  order  of  0.5  to  10 
microns  in  dimensions.  Th^  remain  suspended  in  the  water  for  the 
duration  of  the  absorption  measurements.  The  measured  absorption  of 
water  at  630  nm  as  a  function  of  increasing  concentrations  of  scatterers 
using  models  II  and  V  is  shown  in  figure  14.  Table  1  shows  the 
scattering  coefficient  for  each  quartz  concentration  of  figure  14. 

For  these  measurements  the  calibration  constants  Kx  and  K2  were 
(0.5197  m*^  and  0.8642  m*^)  for  model  n  and  (0.6715  m~^  and 
0.5901  m'^)  for  model  V.  This  measured  absorption  shows  an 
anomalous  Increase  with  increasing  concentration  of  scatterers  beyond 
•1  gm/L  However,  the  effect  is  small;  even  at  2.3  gm/1  the  measured 
absorption  is  onty  in  error  by  -15%.  At  1  gm/1,  where  the  scattering 
has  negligible  effect,  the  extinction  length  due  to  scattering  is  only 
about  7  cm.  Since  the  cavity  diameter  of  model  n  is  3.7  cm.  the 
corresponding  scattering  optical  depth  is  0.53.  At  the  same  time,  due 
to  the  high  reflectivity  of  the  walls  and  the  corresponding  large  number 


of  transits  through  the  sample,  the  effective  extinction  length  due  to 
absorption  by  the  water  is  >-3.7  m. 


Scatterers  (grams/liter) 


X  MODEL  n  o  MODEL  V 

630nm 


Flgnre  14.  Observed  absorption  of  water  at  630  nm  versus 
concentration  of  nonabsorbing  scatterers  in  models  II  and  V. 
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Tdde  1.  Scattering  coefficients  for  various  concentrations  of 
quartz  powder  suspended  in  pure  water. 


Scatterers  (grams/liter) 

Scattering  (1/cm) 
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In  initial  cavity  designs,  the  observed  increase  in  absorption  at 
hi^  concentrations  of  scatterers  was  attributed  to  creation  of  an 
anisotropy  of  the  Fq  irradlance  at  the  sample  surface  by  the  Sq 
detector  fiber.  Specifically,  the  fiber  is  a  radiation  sink  so  the 
radiation  density  in  its  immediate  vicinity  is  decreased.  The  cavity 
design  of  models  n  -  IV  were  intended  to  decrease  this  effect;  the  Fq 
irradlance  is  allowed  to  travel  down  a  Spectralon  tuimel  where  a  fiber 
can  sample  the  intensity  without  directly  contacting  the  sample.  The 
anisotropy  caused  the  detector  will  be  averaged  over  the  larger 
entrance  area  of  this  tunnel.  The  experimental  results  have  shown 
that  this  was  in  fact  not  a  major  source  of  the  problem.  Each  of  the 
design  modifications  of  models  II  -  IV  produced  the  same  response  as 
shown  in  figure  14. 
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Based  on  an  anafytlc  one-dlmenslonal  model  developed  by  Dr.  G. 

W.  Kattawar^®,  we  now  believe  that  the  apparent  increase  in 
absorption  in  the  presence  of  a  hi^  concentration  of  particles  is  due 
to  an  anisotropy  of  the  radiation  density  in  the  sample.  Specifically, 
the  scattering  is  so  great  that  even  with  very  weak  absorption,  the 
center  of  the  sample  is  darker  than  the  region  close  to  the  Spectralon 
wall.  The  scattering  increases  the  effective  path  length  through  the 
sample.  The  effective  path  length  firom  the  sample  boundary  to  the 
sample  midpoint  becomes  long  enough  for  the  absorption  of  the 
medium  (water  in  this  case)  to  appreciably  decrease  the  radiance  in 
the  center  of  the  sample.  Figure  15  shows  the  Irradiance  as  a 
function  of  the  fractional  distance  from  the  sample  boundary  to  the 
center  of  the  sample  for  a  sample  with  an  optical  depth  of  .53  and  a 
cavity  with  dimensions  of  model  n  (based  on  Dr  Kattawar’s  one 
dimensional  model  of  an  ICAM).  This  effect  will  be  dependent  on  the 
concentration  of  scatterers.  the  absorption  of  the  sample,  and  on  the 
diameter  of  the  sample  cell.  If  we  increase  the  cell  diameter  we 
expect  to  see  this  same  effect  at  a  lower  concentration  of  scatterers. 
Figure  14  shows  the  effect  of  increasing  concentrations  of  scatterers 
suspended  in  water  for  model  n  (cell  diameter  of  3.7  cm)  and  for 
model  V  (cell  diameter  of  9  cm).  The  absorption  begins  to  increase  at 
a  lower  concentration  of  scatterers  as  the  cell  diameter  increases. 

This  means  that  at  extremely  high  concentrations  of  scatterers 
our  assumption  of  a  isotropic,  homogeneous  radiance  distribution  in 
the  sample  volume  is  no  longer  valid.  Specifically  the  term  for  the 
energy  absorbed  in  equation  (1)  is  no  longer  valid  at  these  high 


concentrations  of  scatterers  (Fq  Is  not  constant  throughout  the 
sample  volume).  This  does  not,  however,  significantly  limit  the  use  of 
the  ICAM  as  these  concentrations  of  scatterers  are  far  greater  than 
what  would  normally  be  encount^ed.  In  cases  where  it  is  necessary 
to  measure  absorption  in  the  presence  of  such  high  concentrations  of 
scatterers.  an  ICAM  with  a  smaller  sample  volume  should  be  used. 


Fractional  distance  from  boundary 
to  center  of  sample 


Flgnre  15.  Irradiance  vs  distance  from  sample  boundary. 
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B.  Linearity 

linearity  of  the  ICAM  was  tested  by  measuring  the  absorption 
coef^ent  for  a  series  of  (tye  solutions  containing  incremental 
amounts  of  pigment  dissolved  in  filtered/deionized  water.  As 
shown  in  figure  16.  the  data  remained  linear  up  to  an  absorption 
of  25  meter'^  for  model  11.  Models  in  and  IV  were  found  to 
become  nonlinear  at  much  lower  concentrations  of  cfye  than 
models  where  the  quartz  cell  completely  filled  the  diameter  of 
the  inner  cavity.  This  result  is  not  yet  fully  understood. 


Dye  Concentration  (mg^ter) 


Flgnre  16.  Linearity  test;  absorption  versus  dye  concentration. 
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Model  V  remained  linear  up  to  an  absorption  of  15  meter-^.  As 
the  dimensions  of  the  cavity  are  increased  a  nonlinear  result  is 
e:q>ected  at  lower  concentrations  of  dye  due  to  the  longer  path  length 
through  the  absorber.  This  result  indicates  that  accurate  calibration 
constants  can  be  obtained  onty  calibrating  with  samples  which 
avoid  the  nonlinear  region. 

C.  Spectral  Measurements 

Models  II  and  V  have  been  used  to  measure  absorption  spectra  of 
pure  and  natural  waters  and  of  phytoplankton  cultures.  In  August  of 
1989.  model  n  was  used  to  measure  various  samples  at  the  Naval 
Oceanographic  Research  and  Devek>pment  Agency  located  at  Stennis 
Space  Center,  Mississippi.  At  this  time  a  method  for  determining  the 
calibration  constant  K2  without  reference  to  accepted  values  for 
absorption  of  ptire  water  had  not  been  realized.  Since  the  cavity  was 
known  to  measure  total  absorption,  data  firom  Smith  and  Baker^  for 
the  absorption  of  pure  water  was  added  to  the  absorption  of 
calibration  samples  measured  in  a  spectrophotometer.  This  allowed 
us  to  plot  total  absorption  vs  Sj/^  and  to  determine  both  and  K2 
from  this  data  as  described  previously.  This  would  Introduce  an  error 
in  our  measurements  equal  to  the  uncertainty  of  the  Smith  and  Baker 
data  for  absorption  of  pure  water.  This  imcertainty  would  be 
ne^lgible  for  ai^  absorption  significantly  greater  than  that  of  pure 
water  but  would  become  significant  for  very  small  absorptions. 

Figures  17  through  21  reflect  data  taken  after  calibration  using  this 
method.  Figure  17  shows  the  absorption  spectra  of  distilled  water  as 
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well  as  seawater  and  filtered  seawater  taken  firom  the  Gulf  of  Mexico. 
The  difference  between  the  seawater  and  the  seawater  filtrate  spectra 
is  shown  in  figure  18  and  indicates  a  typical  absorption  spectra  for 
marine  phytoplankton  (the  primary  suspended  particulate  found  in 
seawater).  Figure  19  shows  the  difference  between  the  spectra  of  the 
seawater  filtrate  and  distilled  water  and  is  typical  of  a  spectra  for 
dissolved  organic  matter.  Figure  20  shows  the  spectra  for  a  laboratory 
culture  of  phytoplankton  known  as  Dunaliella  tertiolecta.  The  difference 
between  the  spectra  for  the  sample  and  the  filtrate  gives  the 
absorption  of  the  phytoplankton.  To  examine  the  sensitivity  of  the 
ICAM«  a  laboratory  culture  of  Micromonas  pusilla  was  diluted  to  obtain  a 
concentration  of  chloiophyll>a  of  onfy  0.44  mg  per  cubic  meter 
(within  the  range  of  the  purest  natural  waters).  The  difference  in 
absorption  between  the  sample  and  the  filtrate  was  on  the  order  of 
only  .001  meters*^,  yet  the  ICAM  could  resolve  the  characteristic 
spectra  for  this  phytoplankton  species  as  shown  in  figure  21.  The 
fact  that  this  spectra  compares  favorabty  to  the  spectra  for  the  same 
species  measured  in  a  spectrophotometer  using  other  methods 
indicates  that  calibration  errors  due  to  the  uncertainty  of  the  Smith 
and  Baker  data  for  pure  water  is  ne^gible.  Between  400  to  490  nm 
the  spectrophotometer  data  is  abnormal^  high.  The 
spectrophotometer  data  also  shows  an  error  in  the  shape  of  the  curve 
between  4(X)  and  420  nm.  The  absorption  maximum  should  be  at 
420  nm  not  at  400  nm.  This  is  attributed  to  some  procedural  or 
systematic  error  in  the  spectrophotometer  datai^. 
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X  Seawater  •  double  distilled  water 

+  Filtrate 

Flgnre  17.  Absorption  spectra  for  seawater,  seawater  filtrate, 
and  distilled  water. 
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Fignre  18.  Absorption  spectra  for  seawater  particulates. 


Absoi 


400.00  500.00  600.00  700.( 


Wavelength 

X  Filtrate  minus  distilled 


Flgiire  19.  Difference  between  seawater  filtrate  and  distilled 
water. 
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Absorption  spectra  for  DunalieUa  tertiolecta. 
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Model  V  of  the  ICAM  was  completed  in  December  of  1989. 
Calibration  was  accomplished  using  the  method  \(diere  K2  is 
detennined  by  setting  the  absorption  equal  to  zero  with  the  cavity 
empty.  Figures  22  and  23  show  a  comparison  of  several  pure  water 
samples  with  the  data  from  Smith  and  Baker.  Three  samples  of 
filtered-deionized  water  and  one  sample  each  of  HPLC  grade  water 
and  seawater  were  measured  at  various  times  over  a  fotor  day  period. 

The  spectra  agree  well  in  the  wavelength  region  from  550  to  700 
nm  with  the  exception  of  the  points  at  660,  680,  690,  and  700.  This 
discrepancy  is  attributed  to  the  following  calibration  error.  The 
absorption  spectra  for  calibration  sample  1  (shown  in  figure  1 1) 
appears  discontinuous  at  660  and  680  nm  and  the  absorptions  of  all 
three  samples  vary  by  less  than  .006  meters'l  from  690  to  700  nm 
which  would  allow  greater  error  in  determining  Ki  in  this  region. 

In  the  region  frxim  460  to  550  nm  the  ICAM  data  were  consistent 
in  that  the  spectra  for  three  different  samples  of  deionized-filtered 
water  and  HPLC  grade  water  measured  at  different  times  over  a  4  day 
period  all  agree  to  within  an  absorption  value  of  .01  meters'^.  The 
absorption  of  the  seawater  sample  is  consistent  with  the  other 
samples  in  this  region  in  that  the  values  agree  at  the  higher 
wavelengths  and  the  absorption  of  the  seawater  sample  slowly 
increases  as  the  wavelength  decreases,  indicating  some  dissolved 
oxganics  in  the  seawater  sample.  Even  though  the  ICAM  data  is 
consistent  in  this  region  it  is  less  than  the  data  of  Smith  and  Baker 
and  reaches  negative  absorption  values  from  460  to  490  nm  which  we 
know  to  be  erroneous.  This  error  has  not  yet  been  identified  but  may 
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Flgnre  22.  Absorption  spectra  for  various  pure  waters 
compared  to  data  from  Smith  and  Baker.  550  to  700  nm. 
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Figure  23.  Absorption  spectra  for  various  pure  waters 
compared  to  data  from  Smith  and  Baker.  400  to  550  nm. 
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be  due  to  several  causes.  It  may  be  attributed  to  calibration  errors  due 
the  change  in  the  refractive  index  between  the  cavity  empty  and  filed 
with  water.  Specifically,  determining  the  calibration  constant  K2  by 
setting  the  absorption  equal  to  zero  with  the  cavity  empty  assumes  no 
appreciable  change  in  the  optical  properties  of  the  cavity  between  the 
cavity  when  filled  with  air  and  water.  This  error  may  also  be 
attributed  to  contamination  of  the  sample  cell.  If  the  Irgalan  Black 
solutions  used  for  calibration  leave  a  sli^t  residue  on  the  quartz  cell, 
then  the  measurement  of  Sx/Sq  with  the  cell  empty  would  be  higher 
than  if  this  contamination  were  not  present.  This  would  not  affect 
oiir  measurements  as  long  as  the  contamination  remained  constant. 

If.  however,  the  contamination  is  subsequently  removed  due  to 
repeated  rinsing  of  the  sample  cell,  then  all  values  for  Sx/Sq  and 
absorption  would  decrease  by  a  value  equal  to  the  absorption  of  the 
contamination  that  was  removed.  This  implies  a  change  in  the 
intercept,  K2.  of  equation  (3).  If  this  change  is  extremely  small  it 
woiald  onty  become  significant  when  measuring  the  smallest 
absorption  coefficients.  Since  the  absorption  of  pure  water  drops  by 
an  order  of  magnitude  as  the  wavelength  falls  below  550  nm.  it  is 
conceivable  that  the  effect  would  show  up  here.  Furthermore,  if  some 
contamination  were  being  removed  each  time  the  sample  cell  was 
rinsed,  then  each  successive  measurement  would  indicate  a  lower 
absorption  as  more  and  more  of  the  contamination  was  removed.  The 
number  following  each  sample  Identification  in  the  legend  of  figure 
23  indicates  the  sequence  in  which  the  samples  were  measured. 
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Since  successive  samples  of  the  filtered  deionized  water  show 
decreasing  absorption  values,  this  is  a  very  likely  possibility. 

In  the  region  from  400  to  450  nm  the  absorption  values  of  similar 
samples  varied  by  as  much  as  .02  meters*^,  and  were  also  negative. 

The  negative  values  are  attributed  to  the  same  error  as  discussed  for 
the  region  from  460  to  500  nm.  The  variance  in  the  values  is 
attributed  to  noise  of  the  lock-ln  amplifier.  Specifically,  the  gain  for 
the  Sq  signal  in  this  region  is  4  to  10  times  that  used  for  the  rest  of 
the  spectrum.  The  absorption  values  in  this  region  are  also  an  order 
of  magnitude  smaller  than  those  in  the  rest  of  the  spectrum.  This 
leads  to  comparative^  more  noise  in  the  absorption  values,  and 
indicates  that  data  should  be  averaged  over  a  longer  time  period  when 
operating  in  the  lower  wavelength  region. 

C.  Energy  Density  Within  the  Sample 

For  many  absorption  measurements  it  is  desirable  to  know  the 
intensity  of  the  radiation  striking  the  sample.  One  example  of  this  is 
in  measuring  the  absorption  of  phytoplankton.  Absorption  can  be 
affected  if  the  incident  radiation  photoinhibits  the  phytoplankton. 

We  denote  the  scalar  radiance  at  a  point  within  the  cavity  sample 
volume  as  Uf.D)  and  assume  that  it  is  uniform  and  isotropic. 


L(r.D)  =  Lq = constant 


watts 


m^-steradlans  ' 

In  the  cavity  the  radiant  flux  incident  on  the  detector  fiber  from 

within  a  solid  angle  AD  will  be. 

P»L,AQA, 
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where  A  is  the  fiber  area.  The  portion  of  this  power  that  enters  the 
fiber  is  a  function  of  the  angle  of  incidence  at  the  fiber  end.  We  will 
denote  the  fiber  acceptance  function  as  g(6)  where  6  is  the  angle 
between  the  incident  radiation  and  the  fiber  axis.  The  power  in  the 
fiber  will  then  be  Pg(6)  giving  a  signal 

s.=K/pg(e). 

-KL.Ajg(0)dQ, 

a 

=  KL^xJ  gO)  •  sln(0)  •  d0 ,  (5) 

where  K  is  a  constant  which  converts  the  optical  signal  into  units  of 
volts.  We  want  to  find  Lq.  and  solving  equation  (5)  for  it  we  obtain 

L-  ,  So 
®  2icAjK  g(0)  sln(0)  d0  ’ 

We  experimentally  determine  Kg(0)  using  the  experimental  apparatus 
shown  in  figure  24.  The  So  detector  fiber,  coupled  to  its 
photomultiplier  tube,  was  removed  from  the  cavity  and  mounted  on  a 
rotational  platform  with  the  detector  end  of  the  fiber  positioned  at 
the  center  of  rotation.  A  HeNe  laser  of  known  beam  power  was  used 
as  a  source.  The  laser  beam  was  directed  at  the  So  fiber  and  passed 
through  a  neutral  density  filter  and  then  through  a  light  chopper  and 
focusing  lens.  Values  for  So  were  recorded  every  2.5  degrees  as  the 
So  detector  fiber  was  rotated  through  180  degrees  with  respect  to 
the  incident  laser  beam.  This  produced  values  for  So  as  a  function  of 
the  angle  of  incidence  of  the  li^t  source  which  we  call  Sp(0).  We 
assume  the  acceptance  of  the  fiber  is  rotationally  symmetric  about  the 
fiOber  axis.  The  power  incident  on  the  fiber  ^s  Pl  .  the  laser  power 
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after  passing  through  the  neutral  density  filter  and  the  light  chopper. 
The  portion  that  enters  the  fiber  is  P^  gO)  ixdiich  gives  a  signal 

S,(0)*KP,.g(0), 


or 


K.g(0)= 


Sp(0) 
Pl  • 


Substituting  this  into  the  equation  for  Lo  gives 


u- 

4jcAj^  Sp(0)‘  sin(0)  •  d0  , 
we  numerically  integrate  the  Sp(0)  data. 

£^Sy(0)sin(0)d0  =  C. 


to  find 


•  4itAC  * 


The  irradiance  normal  to  any  surface  in  the  cavity  is 


Fo=icLo 


For  a  wavelength  of  628  nm  and  using  the  detector  fiber  from  model 
V,  the  value  for  Lq  is  *048  watts/m^  *  sr.  The  energy  density  (U)  is  Just 

U»4RLQm/c 

where  m  is  the  index  of  refraction  and  c  is  the  speed  of  light. 
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CHAPTER  Vn 

CONCLUSION 

The  ICAM  has  been  shown  to  be  capable  of  measuring  absorption 
independent  of  scattering.  It  can  be  used  to  acciorately  measure  very 
small  absorptions  dess  than  .005  meters*^)  and  to  measure  the 
absorption  of  suspended  particulate  matter. 

There  are  still  some  residual  problems  such  as  identification  and 
correction  of  the  error  that  produced  negative  absorption  values  for 
pure  water  in  the  lower  wavelength  region,  vdiich  will  be  studied  in 
future  work. 
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